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Measurements of double helicity asymmetries in inclusive hadron production in polarized p + p 
collisions are sensitive to helicity-dependent parton distribution functions, in particular to the gluon 
helicity distribution, Ag. This study focuses on the extraction of the double helicity asymmetry in 
r\ production (p + p — > rj + X), the n cross section and the n/n° cross section ratio. The cross section 
and ratio measurements provide essential input for the extraction of fragmentation functions that 
are needed to access the helicity-dependent parton distribution functions. 

PACS numbers: 13.85. Ni,13.88.+e,14.20.Dh 



I. INTRODUCTION 



Until recently, the knowledge about helicity-dependent 
parton distribution functions (PDFs) in the nucleon 
mainly came from next-to-leading order (NLO) QCD fits 
(see, e.g., to the helicity-dependent structure func- 
tion <7i, as measured in fixed-target polarized inclusive 
deep-inelastic scattering (DIS) experiments (see, e.g., 
@> HI)- The resulting helicity-dependent PDF for the 
gluon has rather large uncertainties due to the fact that 
the exchanged virtual photon does not couple directly, 
i.e., at leading order, to the gluon, and that an indirect 
way of accessing it via NLO fits to g\ suffers from the 
limited kinematic reach of the fixed target experiments. 
Accessing the helicity-dependent gluon PDF via the so- 
called photon-gluon-fusion process in semi-inclusive DIS 
has not yet resulted in better constraints, see Ref. 0,i 
and references therein for details. Thus, additional data 
from polarized p + p scattering, in which longitudinally 
polarized gluons are directly probed via scattering off 
longitudinally polarized gluons or quarks, should reduce 
the uncertainties in the helicity-dependent gluon PDF. 
This has been demonstrated in a global NLO fit [6| us- 
ing, for the first time, the available inclusive and semi- 
inclusive polarized DIS data together with first results 
from polarized p + p scattering at the Relativistic Heavy 
Ion Collider (RHIC). The results included were the dou- 
ble helicity asymmetries in inclusive 7r° and jet [T3| 
production from the PHENIX and STAR experiments, 
respectively. 

The double helicity asymmetry in inclusive hadron pro- 
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duction is given as 



J2 abc A/a ® A/fc ® Aa ab ^ cX ' ® D h c 
2a 

where the cross section <j ++ (er H ) describes the reac- 
tion where both protons have the same (opposite) he- 
licity. The helicity-independent cross section is defined 
as a = ((T ++ + a 4, )/2. The helicity-dependent de- 
composition of the numerator is given on the right-hand 
side of Eq. [TJ where A/ Q , A/& represent the helicity- 
dependent PDFs for quarks or gluons, and Act are the 
helicity-dependent hard scattering cross sections calcu- 
lable in perturbative QCD (pQCD). The kinematic de- 
pendences of these terms are omitted for simplicity. At 
leading order (LO) the fragmentation functions can 
be interpreted as the probability for a certain parton c to 
fragment into a certain hadron h and thus they are not 
needed in the case of jet and direct photon production. 
In current global fits of parton helicity distributions, the 
fragmentation functions are assumed to be spin indepen- 
dent. 

This study focuses on the midrapidity cross section 
and double helicity asymmetry in inclusive 77 production 
(p+p — > r\ + X) as a function of transverse momentum 
(pr) and the rj/n cross section ratio at y/s — 200 GeV 
measured at the PHENIX experiment at RHIC. The mea- 
surement of the 77 double helicity asymmetry adds inde- 
pendent data with different systematics to the present 
set of polarized data available to PDF fits. Even when 
compared to a closely related data set, e.g., the PHENIX 
tt° data on double helicity asymmetries, the difference in 
the fragmentation functions can lead to a different sen- 
sitivity to certain helicity-dependent PDFs. In contrast 
to the 7T°, the experimental data available to extract r\ 
fragmentation functions is rather limited. The existing 
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r\ cross section measurements from e + + e~ collider data 
can constrain the quark fragmentation functions to some 
degree, but the extraction of gluon fragmentation func- 
tions requires cither rather precise e + + e~ data taken 
in a wide range of center of mass energies, or cross sec- 
tion measurements from processes where gluons are di- 
rectly involved, e.g., p+p scattering. Therefore, the data 
on cross sections and cross section ratios presented here 
serve as important input for the extraction of fragmen- 
tation functions, in particular as the measurement has 
been performed over a wide range of px- 



II. EXTRACTION OF 77 AND n° YIELDS 

The 77 (tt°) meson is reconstructed through its main 
decay channel, 77 (ir°) — > 77, with a branching ratio (BR) 
of about 39% (99%) [n|. The data were taken at the 
PHENIX [H[ experiment in 2005 and 2006. After data 
quality and vertex cuts, 2.5 pb -1 from the 2005 data and 
6.5 pb -1 from the 2006 data are used for the analysis. 
The data sets from both years are used for the extrac- 
tion of the statistics-limited double helicity asymmetry 
in 77 production while only the larger data set from 2006 is 
used for the extraction of the predominantly systematics- 
limited cross section measurements. Note that the anal- 
ysis described in the following, including all cuts, is done 
in the exact same way for the 77 and the 7r° meson in or- 
der to minimize the systematic uncertainties on the 77/ 7r° 
cross section ratio presented below. 

Two data sets have been analyzed, collected by requir- 
ing two different trigger selections. The minimum bias 
(MB) trigger requires coincident signals in two beam- 
beam counters (BBCs) [l3j], which are arrays of quartz- 
radiator Cerenkov counters providing full azimuthal cov- 
erage at pseudorapidities of 3.0 < |?7| < 3.9. Based on 
the timing of the signals from the two BBCs, the event 
vertex is reconstructed and required to be within 30 cm 
of the nominal interaction point. In addition to the MB 
trigger, the high-p^ triggered data set requires an en- 
ergy deposition larger than approximately 1.4 GeV in an 
area of 4 x 4 towers in the electromagnetic calorimeter 
(EMCal) 0. 

The EMCal is the primary detector used in this anal- 
ysis, located at a radial distance of about 5 m from the 
beam pipe. It covers the pseudo-rapidity range 1 77 1 < 0.35 
and has an azimuthal acceptance of A<p = tt. The EMCal 
consists of eight sectors, six of which are composed of a 
total of 15552 lead-scintillator (PbSc) sandwich towers 
(5.5 cm x 5.5 cm x 37.5 cm), and two sectors of lead- 
glass (PbGl) Cerenkov calorimeters, consisting of a total 
of 9216 towers (4 cm x 4 cm x 40 cm). For the cross 
section measurements only the PbSc was used. 

A cluster in the EMCal is assumed to originate from 
a photon if the following criteria are met. First, since 
showers in the EMCal are not confined to a single tower, 
a shower profile analysis can be used to reject hadrons, 
which usually produce broader showers than photons. 



Since hadrons are slower than photons, an additional 
time of flight cut is used for photon identification. Fur- 
thermore, the cluster must not be associated with a hit 
from a charged particle in the Pad Chamber (PC3) just 
in front of the EMCal; an exception is made if the hit 
position in the EMCal and in the PC3 are aligned in 
such a way that the particle could have come from the 
vertex on a straight line, i.e., it was not bent in the cen- 
tral magnetic field. In this case, the cluster is accepted 
as a photon candidate since it is likely that the original 
photon converted into an e + e~ pair before the PC3 but 
outside the magnetic field. The latter two selection cuts 
are used in the analysis of the double helicity asymmetry 
but not in the extraction of the cross sections, leading to 
a smaller signal to background ratio in the cross section 
measurements. In order to exclude clusters with poten- 
tially incorrectly reconstructed energies due to leakage 
effects, the tower with the largest energy deposition in 
a cluster must not be in the outermost two columns or 
rows of an EMCal sector. In addition, there must not be 
a noisy or dead tower in the eight towers surrounding the 
central one. 

Using all possible pairs of photon candidates, the two- 
photon invariant mass spectrum is calculated. An up- 
per limit of 0.7 is placed on the energy asymmetry, 
\E\ — E 2 \/ [Ei + Ei), of the two cluster energies, E\ and 
E2 , in order to reduce the combinatorial background due 
to numerous low-energy background clusters. It is also 
checked that either of the two clusters coincides with an 
area in the EMCal that caused a high-pr trigger. Fi- 
nally, the pt of the diphoton is required to be larger 
than 2 GeV/c. At smaller pt, the uncertainties in the 
cross section extraction become too large due to large 
backgrounds and limited acceptance. 

The 77 and it cross sections and the 77 double helicity 
asymmetry are extracted in bins of px ■ Using the selec- 
tion cuts for the cross section extraction on the high-px 
triggered data set, the resulting invariant mass distribu- 
tions in the vicinity of the 77 peak are shown in Fig. [1] for 
three different bins of px- For bins at small px, the signal 
extraction is based on fits to the invariant mass distribu- 
tions using a Gaussian for the signal plus a second-order 
polynomial for the background that describe the vicinity 
of the 77 and ir° peaks very well. For px 10 GeV/c 
the signal extraction based on fits becomes unreliable as 
limited statistics leads to large fluctuations in the fit re- 
sults for the mean and width of the peaks. Therefore, 
the mean and width are taken from a Monte-Carlo simu- 
lation, which describes the mean and width of the 77 and 
7T° peaks as a function of px very well for bins below 
px ~ 10 GeV/c, giving confidence in using them for the 
bins above. In these cases, the number of background 
counts under the signal peak is estimated by using the 
number of counts in the sidebands. The sidebands are 
on both sides of the mean of the peak, between 4 and 
7 (4 and 6) times the Gaussian width of the peak for 
the cross section (double helicity asymmetry) analysis. 
However, the exact position and width of the sidebands 



5 



x10 



2.5 GeV/c < p T < 3 GeV/c 




0.7 0.8 
M YT (GeV/c 2 ) 



6.5 GeV/c < p T < 7 GeV/c 




0.8 

M-yy (GeV/C 2 ) 



13GeV/c<p T <14GeV/c 



3 
O 

o 



30 r 
25 r 
20 r 
15 r 
10 - 
5 r 
- 



0.3 



0.4 



0.5 



0.6 



0.7 0.8 

M Y7 (GeV/c 2 ) 



FIG. 1: Invariant mass distributions in the vicinity of the r\ 
peak for the high-pT triggered data set and for three different 
bins of pt- The selection cuts for the cross section extraction 
are used. 



are varied and possible effects are taken into account in 
the systematic uncertainty. In the mid-pr range between 
3 GeV/c and 10 GeV/c, where statistics is sufficient for 
the fit results to be stable and the background in the 
vicinity of the peaks is approximately linear so that the 
sideband subtraction is applicable, both methods agree 
as expected. 



III. THE r) CROSS SECTION AND V /ty° CROSS 
SECTION RATIO 

The r\ and 7r° meson cross sections are calculated from 



E 



,d 3 cr _ 111 1 N(Apr, Ay) 
d 3 p 2npT BR C A e t rig e rec Apr Ay 



(2) 



where C denotes the integrated luminosity, A the accep- 
tance, etrig the trigger efficiency, e rcc the reconstruction 
efficiency, and N the number of reconstructed mesons. 

The luminosity is calculated from the number of MB 
events divided by the cross section for events selected by 



the MB trigger. For the latter, a value of 23.0 mb with 
a systematic uncertainty of 9.7% has been derived from 
Vernier scan results [l5| and an extrapolation for subse- 
quent years. The acceptance is calculated from a Monte- 
Carlo simulation using, as an input, the map of noisy and 
dead cells also used in the data analysis. The systematic 
uncertainty on the acceptance calculation is 3.6% (3%) 
for the rj (ir°) meson. The trigger efficiency for the MB 
data is given by the MB trigger efficiency for rj and 7T° 
production. The trigger efficiency for the high-pT trig- 
gered data set is given by the MB trigger efficiency times 
the efficiency of the high-pT trigger. The MB (high-py) 
trigger efficiency is determined by the ratio of the number 
of reconstructed ir° or rj mesons in a high-pT (MB) trig- 
gered sample in coincidence with the MB (high-py) trig- 
ger divided by the number of reconstructed 7r° or rj events 
without the coincidence. The MB trigger efficiency is 
0.78 for both tt° and r\ mesons over the whole range of 
Pt considered, with a systematic uncertainty of 3%. The 
high-pT trigger efficiency reaches a plateau at a level of 
0.90 for pt > 4.5 GeV/c and is also very similar for both 
mesons. Due to the fact that the turn-on curve of the 
high pt trigger is very steep and reaches an efficiency of 
about 0.80 in the 3 < pt < 3.5 GeV/c bin, the efficiency 
for pt < 3 GeV/c has a large systematic uncertainty. 
Therefore, the cross section calculation is based on the 
smaller MB triggered data set for pt < 3 GeV/c and on 
the high-pT triggered data set for larger transverse mo- 
menta. The reconstruction efficiency accounts for loss of 
photons due to conversion (6% ±2%) and due to the cut 
on the shower shape discussed above (4% ±2%). In the 
case of 7T° production, merging of the two decay photons 
into a single cluster is considered for pt > 10 GeV/c. 

The rj cross section as a function of px between 2 and 
20 GeV/c is shown in Fig. [5] and tabulated in Table Q] 
Note that a bin-shift correction is applied in order to be 
able to plot each data point at the center of each given 
Pt bin, which, due to the exponentially falling spectrum, 
does not represent the true physical value of the yield in 
that bin [16[ . This in particular facilitates the calculation 
of the tj/tt cross section ratio. 

The rj cross section is consistent with an earlier 
PHENIX measurement [ijj covering a smaller range in 
Pt from 2.5 to 12 GeV/c. The error bars shown in Fig. [2] 
are the statistical and systematic uncertainties added in 
quadrature. Not included is an overall normalization un- 
certainty of 9.7% due to the uncertainty in the luminosity 
measurement. The other dominant systematic uncertain- 
ties are an approximately p^-independent uncertainty of 
about 8% due to the uncertainty on the global energy 
scale of 1.2%, possible non-linearities in the energy scale 
affecting mainly points with px > 10 GeV/c, and uncer- 
tainties from the signal extraction affecting principally 
the two lowest pt points, which have a large background 
underneath the 77 peak. The systematic uncertainties 
are subdivided into uncertainties that are uncorrelated 
between pt bins (type-^4), correlated between px bins 
(type--B), and overall normalization uncertainties (type- 
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TABLE I: Results, statistical, and systematic (type-_Bi, type- 
B2) uncertainties of the measured rj cross sections from the 
2006 data set. There is an additional normalization uncer- 
tainty of 9.7% (type-C). 



9.7% normalization uncertainty 
is not included 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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FIG. 2: Cross section for midrapidity inclusive 77 production 
at ^fs = 200 GeV as a function of pr and its comparison to 
NLO pQCD calculations at three different scales /j,. The er- 
ror bars shown are the statistical and systematic uncertainties 
added in quadrature. Not included is the overall normaliza- 
tion uncertainty of 9.7%. Note that the fragmentation func- 
tions used in the calculations are partially constrained by this 
data. See text for details. 
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C). As described above, the peak extraction is based on 
different methods depending on px- Thus, the px bins in 
certain regions are correlated, but there is no full correla- 
tion over the whole range. Such kind of uncertainties are 
sub-categorized as type-i?i, in order to distinguish from 
those correlated over all pt bins (type-i^)- All other 
uncertainties, except the one from the luminosity mea- 
surement (type-C), are assumed to be in this category. 

The 77 cross section from p + p scattering presented 
here, together with the above mentioned earlier PHENLX 
data in a smaller range in pr , and various r\ cross section 
measurements from e + + e~ scattering have been used in 
a global fit to extract new fragmentation functions for r\ 
production at NLO [l8j]. Earlier determinations of r\ frag- 
mentation functions based on SU(3) model estimates at 
LO and normalizations taken from a Monte Carlo event 
generator at NLO are described in Refs. [H, [2(1 and 
Ref. " 
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respectively. Due to the absence of data on 
semi-inclusive 77 production the fragmentation functions 
can only be extracted separately for each quark flavor 
with additional assumptions. The assumption that all 



light quark fragmentation functions are the same, i.e., 
Dl = D% = = Dl = D\ = D], has been used 
in Ref. Using these fragmentation functions and 

the CTEQ6M [H PDFs as an input to the NLO code 
of Ref. [23[ , pQCD calculations at three different scales 
\i are carried out. Here, /1 represents the factorization, 
renormalization and fragmentation scales, i.e., the three 
scales are set equal in each separate calculation. With 
these new fragmentation functions, for which the present 
data constitute nearly 20% of the input experimental 
data points, the cross section is described well. 

The contributions of the various scattering subpro- 
cesses, gluon-gluon (gg), quark-gluon (qg), and quark- 
quark (qq), to the 77 production as a function of pr, are 
shown in Fig. [3] For comparison, they are also shown 
in the case of n° production [24|. While the correspond- 
ing uncertainties are difficult to quantify, it is clear that 
the subprocess contributions to the r\ and ir° production 
are, within uncertainties, identical up to a px of approx- 
imately 10 GeV/c. This is the kinematic range of the 77 
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FIG. 3: (color online) Fractional contribution of gluon-gluon 
{gg), quark-gluon (qg), and quark-quark (qq) scattering to 
the r\ production in the pQCD calculation of Fig. [2] and to 
the 7T° production [Sj, as a function of pr- 



and 7T° double helicity asymmetries presented below and 
published in Ref . @, H| , respectively. Consequently, these 
measurements have approximately the same sensitivity to 
the gluon helicity distribution accessible via the gg and 
qg subprocesses. The differences at larger values of px 
are mostly related to uncertainties in the fragmentation 
functions and thus do not necessarily indicate different 
sensitivities to polarized PDFs. 



function of pr, is given in Fig. [4] and Table [Til The ratio 
has been extracted in a single pass over the same data 
set, thus minimizing systematic uncertainties. In partic- 
ular, the large normalization uncertainty of 9.7% arising 
from the luminosity calculation cancels completely. Also 
all other systematic uncertainties are assumed to either 
cancel or be reduced to a negligible amount with the ex- 
ception of the following. The systematic uncertainties 
due to the -q and ir° peak extraction (type- Si) and due 
to the correction for possible merging of the two 7r° de- 
cay photons into a single cluster (type-i?2) do not can- 
cel. Furthermore, while the uncertainties on the high-p^ 
trigger efficiency (type-i?2) are assumed to cancel above 
Pt = 4.5 GeV/c where the efficiency is flat, a remaining 
2% uncertainty on the ratio is assigned for differences 
in the trigger turn-on curve for 3 < px < 4.5GeV/c. 
Finally, the systematic uncertainty on the acceptance 
(type-i^) is reduced to a px independent contribution 
of 2%. The ratio is presented up to px — 14 GeV/c 
only, as beyond this point the statistical and systematic 
uncertainties become rather large. The latter is due to 
the fact that for increasing transverse momenta the two 
photons from the 7T° have a strongly increasing probabil- 
ity of being reconstructed as only a single cluster in the 
calorimeter, leading to a rather large systematic uncer- 
tainty arising from the correction for this effect. 

Except for an initial increase due to the different me- 
son masses, the data do not exhibit a strong dependence 
on pt- This hints towards a similar dependence of the r\ 
and 7T° fragmentation functions on the energy fraction of 
the parton carried by the hadron. Fitting the cross sec- 
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FIG. 4: Cross section ratio for the midrapidity production 
of inclusive r\ to ty° mesons at yfs — 200 GeV as a function 
of pt- The error bars show the statistical and systematic 
uncertainties added in quadrature. The solid curve shows the 
ratio of the NLO pQCD calculations shown in Fig. [3] and 
the corresponding one for the 7r°. The dashed curve shows 
the result of a pythia Monte-Carlo simulation. See text for 
details. 

Constraining the fragmentation function further 
should be possible by including precise rj to ir° cross sec- 
tion ratios in the extraction. The cross section ratio, as a 



TABLE II: Results, statistical, and systematic (type-Bi, 
type-i?2) uncertainties of the measured r\ to tv° cross section 
ratio from the 2006 data set. 
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tion ratio, including its statistical and type-i?i system- 
atic uncertainty, to a constant, gives R^j^o = 0.51 ±0.01 
(x 2 /ndf = 18.3/17), with a remaining systematic uncer- 
tainty of 0.01 from the type-i?2 systematic uncertainty. 
This result does not change when fitting the data above 
p T = 3 GeV/c (x 2 /ndf = 14.9/15) instead of fitting the 
full range. 

Within the uncertainties, the present measurement of 
77/71"° is consistent with all previous measurements in p+p 
collisions, going back to the measurement reported in 
Ref. [25| . A detailed com par ison of subsequent measure- 
ments is summarized in (l7| . The observed ratio is in 
good agreement with a pythia 6.131 [26[ calculation [ijj 
shown in the same figure, which is using the default set- 
tings and the Lund string fragmentation model [27| ■ The 
solid line in Fig. @] shows the ratio of the NLO pQCD 
calculations at a scale \l — pr (see Fig. [2]) and the corre- 
sponding one for the 7T° using the same PDF but the 7r° 
fragmentation function of Ref. [24[ ■ Note that the shape 
of this calculated cross section ratio is not necessarily 
well determined as the statistical uncertainty on the 77 
fragmentation function, defined by Ay 2 = 2%, results in 
an uncertainty on the 77 cross section between about 5% 
and 9%, depending on pt [Hj]. 

The calculated ratio underestimates the data even 
though the 77 cross section presented in this paper and 
earlier PHENIX tt° data are part of the input in the ex- 
traction of the fragmentation functions. This indicates 
that the constraints from the separate fits are less strin- 
gent than fitting the cross section ratio directly. This 
is obvious from the fact that some of the experimental 
systematic uncertainties cancel in the ratio as already 
discussed above, in particular the overall normalization 
uncertainty of 9.7% due to the uncertainty in the lumi- 
nosity measurement. For example, the earlier PHENIX 
data used in the extraction of the 7r° fragmentation func- 
tions was scaled by a factor of 1.09 [24| in the fit, which 
is within the experimental normalization uncertainty, but 
leads to a smaller calculated cross section ratio as can be 
seen in Fig. |U Also, the dependence of the calculated 77 
and 7T° cross sections on the theoretical scale, as shown 
in, e.g., Fig. 2 largely cancels in the calculation of the 
ratio [201 ] . Hence it appears that improved constraints on 
77 and 7r° fragmentation functions can be derived by di- 
rectly including the data on the 77/71"° cross section ratio 
in the fit. 



IV. DOUBLE HELICITY ASYMMETRY FOR 77 
MESONS 



Experimentally, the double helicity asymmetry (Eq. Q} 
translates into 



All = 



1 



RN 4 



Pb Py 



RN 4 



with R = 



(3) 

where N ++ (JV+_) is the experimental yield for the case 
where the beams have the same (opposite) helicity. The 



polarizations of the two colliding beams at RHIC are de- 
noted by Pb and Py- The relative luminosity R is mea- 
sured by a coincident signal in the two BBCs that sat- 
isfies the vertex cut. Uncertainties on All of 2 x 10~ 4 
(7 x 10~ 4 ) for the 2005 (2006) data due to relative lumi- 
nosity uncertainties are uncorrelated between years. The 
asymmetries and uncertainties are combined by weight- 
ing by all year-to-year uncorrelated uncertainties in each 
Pt bin. The results are given in Table IIIII 

The degree of polarization is determined from the com- 
bined information of a polarized-proton on carbon (pC) 
polarimeter |28(, using an unpolarized ultra-thin carbon 
ribbon target, and from elastic p + p scattering, using a 
polarized atomic hydrogen gas-jet target [29}. The av- 
erage polarization value for the data from 2005 (2006) 
is 0.49 (0.57). There is a relative uncertainty of 4.8% 
in the product of the beam polarizations, correlated be- 
tween the 2005 and 2006 data sets, which is a scale un- 
certainty on the combined asymmetry result, affecting 
both the central values and the statistical uncertainties 
such that the statistical significance of the measurement 
from zero is preserved. Uncertainties on the products 
of the beam polarizations that are uncorrelated between 
years are combined using the same weight factors as for 
the uncertainties due to relative luminosity, and given 
in Table IIIII In order to avoid false asymmetries due to 
a possible variation of detector response versus time or 
due to a possible correlation of detector performance with 
the RHIC bunch structure, all four helicity combinations 
in the colliding bunches are present within four consec- 
utive bunch crossings. Possible transverse components 
of the beam polarizations at the PHENIX interaction 
point are monitored by measuring the spin dependence of 
very forward neutron production [3(| in the zero degree 
calorimeters [3L| . The average transverse component of 
the product in the 2005 data set is less than 0.014±0.003, 
described in more detail in Q, and was measured to be 
negligible in the 2006 data set. 

The double helicity asymmetry for 77 production, 



A\l 



A 



V+BG 
LL 



1 -r 



with 



N BG 



NV + N BG ' 



(4) 



can be calculated by separately measuring the asymme- 
try in the 2er window around the mean of the 77 peak 
(A L + L BG ) and in the sidebands {A B< £) defined above. 

The latter is consistent with zero. The resulting back- 
ground corrected asymmetry for 77 production as a func- 
tion of pt from the combined 2005 and 2006 data is shown 
in Fig. [5] and tabulated in Table Hill It is consistent with 
zero over the measured range as can be expected based 
on the similar contributions of the various scattering sub- 
processes to the 77 and 7r° production shown in Fig. [3] and 
the fact that the double helicity asymmetry for 7r° pro- 
duction Q is consistent with zero as well. 

As can be seen in Fig. [5J the 77 double helicity asymme- 
try is in agreement with NLO pQCD calculations using 
the above mentioned fragmentation functions and two 
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FIG. 5: Double helicity asymmetry for midrapidity inclu- 
sive 77 production from the combined 2005 and 2006 data 
at y/s = 200 GeV as a function of pr- The grey boxes 
are point-to-point systematic uncertainties due to polariza- 
tion and relative luminosity uncertainties and are correlated 
point-to-point, moving all points in the same direction but 
not by the same factor. An additional systematic uncertainty 
of 4.8% on the vertical scale due to the uncertainty in the 
beam polarizations is not shown. The results are compared 
to NLO pQCD calculations using two different sets of polar- 
ized PDFstg, 132]. See text for details. 



different sets of polarized PDFs[6, 32j as an input to the 
code of Ref. 

These data can be used in global fits in order to fur- 
ther constrain polarized PDFs, in particular the helicity- 
dependent gluon PDF. In the future, with improved 
statistics and the availability of flavor-separated fragmen- 
tation functions, double-helicity asymmetries in rj pro- 
duction can potentially constrain the polarized strange 
quark PDF (As) due to the additional s-quark contribu- 
tion in the r\ wave function. Special interest in As arises 
from the fact that its value is negative, when extracted 
from analyses of inclusive DIS data, using hyperon de- 
cay data and assuming SU(3) flavor symmetry 0, Hj], 
but consistent with zero, when directly extracted from 
semi- inclusive DIS data [33- 35] . Global fits can con- 
strain PDFs by simultaneously describing a wide variety 
of experimental channels over a range of kinematics with 
different sensitivities, different experimental systematic 
uncertainties, and different sources of theoretical uncer- 
tainty. Thus, the data presented here open up a valuable 
new channel to improve knowledge of polarized PDFs. 



V. SUMMARY 

The double helicity asymmetry in r\ production is mea- 
sured and found to be consistent with zero in the trans- 
verse momentum range between 2 GeV/c and 9 GeV/c. 
The r\ cross section is determined over seven orders of 
magnitude between 2 GeV/c and 20 GeV/c in trans- 



TABLE III: Double-helicity asymmetry values and uncertain- 
ties for combined 2005 and 2006 data sets. Systematic un- 
certainties given are type B2, scaling all points in the same 
direction but not by the same factor, and are due to polariza- 
tion (a^ y3 ) and relative luminosity (o-f ya ) uncertainties that 
are uncorrelated between years. There is an additional type 
C systematic uncertainty of 4.8% on the vertical scale due 
to uncertainty in the beam polarizations that is correlated 
between years. 



( P t) (GeV/c) 


A LL 


0"stat 


p 

^syst 


R 

V syst 


2.33 


-0.0073 


0.0055 


0.0006 


0.0004 


3.35 


0.0000 


0.0070 


0.0000 


0.0004 


4.38 


-0.0026 


0.0124 


0.0002 


0.0004 


5.40 


0.0137 


0.0228 


0.0016 


0.0004 


6.74 


-0.0111 


0.0320 


0.0021 


0.0004 



verse momentum. In particular due to the wide range 
in transverse momentum these data serve as important 
input for the extraction of fragmentation functions. The 
r\ to 7T° cross section ratio as a function of px has been 
extracted in a single pass over the same data set, thus 
minimizing systematic uncertainties. A fit to a constant 
above pt — 2 GeV/c or px = 3 GeV/c yields a value 



of R 



r)/7T Q 



= 0.51 ± 0.01 stat ± 0.01 syst . The inclusion 



of these data on the ratio in future fragmentation func- 
tion extractions should allow for more precise results for 
both particle species. This opens up the possibility to 
include the data on the double helicity asymmetry in fu- 
ture NLO pQCD fits in order to further constrain the 
polarized parton distribution functions, in particular the 
hclicity-dependent gluon distribution function. 
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